Exchange bias driven by the Dzyaloshinskii-Moriya interaction and ferroelectric 
polarization at G-type antiferromagnetic perovskite interfaces 
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Exchange bias is usually rationalized invoking spin pinning effects caused by uncompensated 
antiferromagnetic interfaces. However, for compensated antiferromagnets other extrinsic factors, 
such as interface roughness or spin canting, have to be considered to produce a small uncompensa- 
tion. As an alternative, here we propose two (related) possible mechanisms, driven by the intrinsic 
Dzyaloshinskii-Moriya interaction and ferroelectric polarization, for the explanation of exchange 
bias effects in perovskites with compensated G-type antiferromagnetism. One of the mechanisms is 
only active when a multiferroic material is involved and it is controllable by electric fields. 

PACS numbers: 75.70.Cn, 75.10.Hk, 75.30.Et, 75.80.+q 
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Introduction. The exchange bias (EB) effect, charac- 
terized by a shift of the magnetic hysteresis loops away 
from the center of symmetry at zero magnetic field, is 
widely reported to exist in magnetic systems where there 
is an interface between antiferromagnetic (AFM) and fer- 
romagnetic (FM) (or ferrimagnetic) materials [l|. 

Theoretically, the EB is understood as induced by spin 
pinning effects at the FM/AFM interface. An uncompen- 
sated AFM interface is usually invoked to illustrate how 
the pinning may work. These uncompensated AFM spins 
at the interface are expected to pin the nearest-neighbor 
(NN) FM spins via the exchange coupling, giving rise to 
a preferred direction for the FM moments. However, de- 
spite its physical appeal, this simple picture is not enough 
to fully understand several real EB cases in a variety of 
magnetic systems. This approach usually predicts an EB 
larger than measured, and also fails to answer why there 
is EB in some fully compensated AFM interfaces 0|- 

Precisely for the subtle case of compensated AFM in- 
terfaces, extrinsic factors are also often considered, such 
as interface roughness |3|. Spin-canting near the inter- 
face can also contribute to the EB Other models 
have also been proposed, such as frozen interfacial and 
domain pinning. Most of these models still need a small 
"frozen" uncompensation of the AFM moments near the 
interface, thus remaining under much debate 0- 

Recently, remarkable improvements in oxide thin-film 
techniques have allowed for the growth and character- 
ization of complex oxide heterostructures with (near) 
atomic precision, opening an avenue for the fabrica- 
tion of multifunctional devices using strongly correlated 
electronic materials [5|. In this context, EB has been 
observed in BiFeOs (BFO) based heterostructures 



More interestingly, the EB in multiferroic heterostruc- 
tures is widely believed to be controllable by electric 
fields. In addition, the EB has also been observed in 
SrRuOa/SrMnOs (SRO/SMO) superlattices 0. Consid- 
ering that both BFO and SMO are well-known compen- 
sated G-type AFM materials (all NN spins are antiparal- 
lel) and that the interfaces are very smooth, the origin of 
the EB in these heterostructures remains a puzzle. The 
purely magnetic interactions framework stemming from 
traditional metallic magnetism appears incomplete to 
deal with the complex physics unveiled in these strongly 
correlated electronic systems, and to address the prac- 
tical matter of how to control the EB by electric fields 
when a multiferroic material is involved. Therefore, new 
mechanisms that emphasize the many simultaneously ac- 
tive degrees of freedom in correlated electron systems are 
needed to better understand these interesting effects. 

The model. Here, we propose two (related) mecha- 
nisms for EB generation in interfaces involving FM/G- 
AFM perovskites. In these mechanisms, the G-AFM in- 
terface can be fully compensated, namely the tiny un- 
compensation caused by various uncertain factors is no 
longer essential (although it can still exist). Therefore, 
our proposed mechanism is conceptually different from 
ideas based on tiny frozen uncompensated AFM moments 
@, I3|- Instead, the interactions between spins and lattice 
distortions become the key intrinsic driving force for the 
mechanisms presented below. 

Let us start with the spin-spin interaction in per- 
ovskites, with a Hamiltonian: 
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where Jij is the standard superexchange (SE) couphng 
between NN spins; i and j are site indices; and S are 
spin vectors. For several large-spin transition metal 
cations in perovskites, such as Mn'^+ and Fe^+, adopt- 
ing the widely-used classical approximation is reason- 
able. In the following, the normalization |5| = 1 will 
be used (the actual magnitude S of the spins can be ab- 
sorbed in a redefinition of couplings) . The second term is 
the Dzyaloshinskii-Moriya (DM) interaction, which arises 
from the spin-orbit coupling 0]. Since \D\ is much 
smaller (by two or three orders of magnitude) than J Q , 
the DM interaction is often neglected. Originally, the 
DM interaction was introduced to explain the presence 
of weak ferromagnetism in AFM materials because the 
DM term can produce a small spin-canting. Recently, the 
DM interaction has also been highlighted as the origin of 
a finite ferroelectric (FE) polarization (P) in multiferroic 
materials with spiral spin order [10]. 

In perovskites, the DM interaction is determined by 
the oxygen octahedron tilting. Usually, the A-site cations 
in perovskites are too small to maintain a stable cubic lat- 
tice. Then, the oxygen octahedra surrounding the B-site 
cations will tilt for a closer packing ll|. The tilting can 
be characterized by the Glazer notation: e.g. a~b~c^ 
where the three letters denote the rotation angle ampli- 
tudes about the [100], [010], and [001] axes, respectively; 
the positive (negative) superscript indicates that the ro- 
tations of two neighboring octahedra, along the tilting 
axis, are in the same (opposite) direction [12l ]. For in- 
stance, in the orthorhombic lattices (e.g. bulk LaMnOa 
at low temperature (T)), the tilting a~a~6^ receives the 
name "GdFeOs-type distortion" and it corresponds to 
rotations around the [110] (dominant) and [001] (sub- 
dominant) axes of the cubic unit cell [l3|. For the M-0- 
M bond (M: B-site metal and O: oxygen), this octahe- 
dral tilting moves the oxygen anion perpendicularly away 
from the midpoint between NN metal cations, as shown 
in Fig. 1(a). Since the tilting rotation is collective, the 
NN oxygens in the same direction (Oi and O2 in M-Oi- 
M-O2-M) should move in opposite directions, namely the 
NN displacements are staggered. 

DM driven EB. From symmetry argumentations, the 
Di^j vector should be perpendicular to the M^-O-Mj bond 
0, as shown in Fig. 1(a). Thus, the D vectors between 
NN bonds along the same direction are also staggered, 
namely Di^i^i——Di^i^i^2- To simplify the discussion, 
let us consider the case where the rotations of NN oc- 
tahedra along the [100] and [010] axes are in opposite 
directions, namely a~b~c* (* can be — , or 0). 

For simplicity, all spins in the AFM and FM side are 
assumed to be collinear. However, because of the differ- 
ent easy magnetic axes or planes for different materials, 
in general the NN spins are noncoUinear at the FM/ G- 
AFM interface (Fig. 1(b)). There are two vectors that 
are staggered: (1) the AFM interface spins Sf-^^ given 
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FIG. 1: (Color online) (a) The (mutually perpendicular) re- 
lationship between the Mi-O-Mj bond, oxygen displacement, 
and Di.j vector, (b) Sketch of the interface between FM and 
G-AFM perovskites, including the oxygen octahedral tilting. 
The staggered directions of the Dij vectors at the interface 
are marked as in- and out-arrows, while the uniform ho vec- 
tors are also shown near the oxygens, (c) The uniform ho 
should be perpendicular to S^^^ and D. 

by (— l)*^^, and (2) the Dij vectors across the interface 
given by {—1/3, where i denotes the site (or bond) se- 
quence at the (001) interface. Combining these two stag- 
gered components Dij and Sf^^^ , it is straightforward to 
obtain a uniform DM effect at the interface: 
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where denotes the spin at the FM side and i and j 
only sums over the interface, /id is the effective magnetic 
field that points into the direction D x (Figs. l(b,c)). 
Note that /id is uniform and independent of the FM 
spins' direction. The combination of Di j and Sf^^, 

namely /id, can be fixed by the field-cooling process and 
then assumed to remain frozen at low T during the hys- 
teresis loop measurement [l3]. Thus, this provides a bias 
field caused by the DM interaction which can produce a 
EB at interfaces of FM/G-AFM perovskites. 

FE driven EB. In the previous discussion, the second 
term (DM interaction) of Eq.(l) was proposed as the mi- 
croscopic origin of EB in generic FM/G-AFM oxide het- 
erostructures. However, the first term (SE) can also con- 
tribute to the EB if multiferroic materials are involved 
in the heterostructure. In ferroelectric (FE) materials, 
spontaneous relative displacements between cations and 
anions induce an electric polarization. Consider the oxy- 
gen positions at the interface shown in Fig. 2(a): in ad- 
dition to the previously mentioned staggered displace- 
ments, that do not induce a finite P, in some multi- 
ferroic materials their FE properties can be assumed to 
be caused by additional displacements of NN oxygens 
that should all be along the same direction to avoid a 
global cancellation. Therefore, the bond-angles at the in- 
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FIG. 2: (Color online) (a) FE-polarization-driven asymmet- 
ric bond-angles and modulated normal SE at the interface. 
A switch of the FE polarization from left to right will also 
switch hj. (b) The estimated |hj| as a function of XpE for 
different values of Xo (from to 0.03 A). The lattice constant 
is assumed to be 4 A and Jj^ is in proportional to cos*(0i,j) 
where 6 is the bond angle All displacements (in units of 
A) are assumed to be coplanar for simplicity. 

terface can become asymmetric by the simultaneous con- 
sideration of these two displacement modes. Since the 
SE coupling magnitude is dependent on the bond-angle 
[l^ . these modulated bond-angles induce an interfacial 
SE coupling J that is also staggered, with values that 
are denoted here as Jl and Js- Once again, as with 
the DM-driven EB, two staggered effects (alternating SE 
couplings at the interface, and alternating spin orienta- 
tions on the AFM side of the interface) compensate each 
other. By this procedure, it is straightforward to obtain 
an additional uniform effective field at the interface: 



interface 
SE 
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Here, /ij — —SjS^ is the effective magnetic field, where 

= {Jl — •^s)/2- When an electric field is applied par- 
allel to the interface to change the uniform polarization 
P, the hj will change simultaneously, namely it is an 
electric- field- controllable EB which is potentially impor- 
tant to design multiferroics devices. 

Both hj and /id may have components parallel to the 
measuring field, although they are perpendicular to each 
other. Experimentally, by varying the electric-field direc- 
tion, estimations for the components of hj and Ed can 
be obtained separately, since Hd is almost independent 
of the FE P, in a first-order approximation . 

Discussion. The basic physical picture related to the 
proposed DM- and FE-driven EB appears clear, but there 
are several practical issues that should be addressed. 

First, in the derivations above, both mechanisms are 
independent of the details of the FM spins. Therefore, 
both mechanisms should be valid for a variety of FM ma- 
terials such as perovskites or metallic alloys @. The 
only condition needed is that the oxygen octahedra of the 
interfacial AFM cations must be complete, i.e. oxygens 
must bridge the two materials at the interface. 

Also, our model should be robust against other tilting 
modes. For a general a'^b^c* mode, the NN -Dy 's at the 



interface are not uniform as long as a and (3 are not 
both simultaneously zero. If this is the case, a net ho 
is still induced, with direction and value varying with 
the mode. For the tilting mode which only rotates along 
the [001] axis (the a^a^c'^ mode in the perfect tetragonal 
lattice), the DM contribution at the (001) interface is 
zero. In this case, /ij will also be zero since the bond- 
angles are uniform. However, there is evidence that many 
perovskite films are not perfectly tetragonal (iBj . 

In addition, since the DM coupling is very weak (par- 
ticularly in nearly tetragonal thin films), it is necessary 
to check whether the EB that it generates is compati- 
ble in magnitude with the experimental EBs. Consid- 
ering ft,D to be only effective at the interface while the 
external magnetic field is applied on all FM spins, the 
maximum EB (when the measuring field is collinear with 
/id) can be estimated as: h-E.B~\hT)\l d—H^^'^^'^'^'' / {dm), 
where d is the FM material thickness in unit cells, and 
TO is the magnetic moment of the FM cation. In a first- 
order approximation, \Di_j\ is proportional to the oxy- 
gen displacement Xq: H^l^'^^'^''^-fXo, with_y the DM 
coefficient roughly estimated as 1 meV/A [lO|. A tiny 
distortion of the M1-O-M2 bond across the interface, as 
small as a 1° bend [i3|, can result in i7^'j^"f'^'='=«0.0175 
meV if the lattice constant is 4 A, indeed very weak com- 
pared with J which is usually larger than 10 meV for per- 
ovskites. Assuming typical values d=lO and to=3 Bohr 
magnetons, the DM driven EB is 100 Oe which is of the 
same order of magnitude as the experimentally measured 
EBs in perovskite heterostructures [7] . 

For perovskite heterostructures involving multiferroics, 
both /iD and hj should be considered. The estimated 
|/ij| vs. the FE oxygen displacement (Xfe, which is 
proportional to the in-plane projection of P) and para- 
metric with Xo are shown in Fig. 2(b). When both 
Xq and X-pE are small, \hj\ behaves approximately as 
XJXqXfe (Aw3.7-4.0). Thus, \hj\/\hu\ is estimated to 
be AJXfe/7, which may be larger than 1 in BFO. 

It is also important to analyze if the DM- and FE- 
driven EB mechanisms are robust against roughness, 
which often is appreciable at interfaces, although recent 
experimental progress in thin-films substantially reduces 
this extrinsic effect. Since there are several uncertain fac- 
tors controlling the interfacial roughness, it is difficult to 
reach robust conclusions from the theoretical perspective. 
For this reason a simplified analysis will be given here, by 
assuming that the FM and AFM cations can be mixed 
near the interface but they not diffuse into inner regions, 
as shown in Fig. 3. If the G-AFM spin order is stable 
enough and there are no crystal defects at the interface, 
the D vectors across the interface alway change simul- 
taneously with the corresponding spins in the G-AFM 
side following the roughness. In other words, the rough- 
ness geometry would not change the combination of Dij 
and G-AFM spin vectors at the interface. Even though 
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FIG. 3: (Color online) Sketch of atomic-scale interface rough- 
ness. Only the ideal G-AFM spin order is shown by arrows. 
The alternation of the Dij vectors across the interface are 
shown as in- and out-arrows. In addition, in the roughened 
case, the D vectors of the (100) and (010) bonds (open/full 
squares) will also be active for the EB. 

the AFM spins can be canted at the roughened regions, 
this may decrease but will not cancel the global /id, as 
long as there are no separated 180° magnetic domains 
or ferroelastic walls. Similarly, it can be shown that the 
FE driven mechanism will not be canceled by roughness 
either. Therefore, both the DM- and FE-driven mecha- 
nisms for EB should in principle work, even in the pres- 
ence of weak interface roughness. 

Note also that the DM- and FE-driven EB are 
anisotropic (related to the crystal direction). Ideally, if 
the measuring field is applied perpendicular to /id (if no 
multiferroics are involved), there would be no EB. A pos- 
sible example is the case of LaMnOa/SMO superlattices, 
in which no EB has been observed using an in-plane mea- 
suring field since all spins are almost in-plane coUinear 
(thus hi) is out-of-plane) [l3|. 

Finally, we remark that we have tested our argumen- 
tations using numerical techniques on a heterostructure 
[lij , and a robust EB in the hysteresis loop was obtained 
by considering our two mechanisms. In the simulation, 
spin canting effects (that can originate from finite-T fluc- 
tuations, exchange couplings at the interface, or magnetic 
field reorientation) are included, but they are not found 
to affect our results qualitatively. 

Conclusions. Here it was proposed that both the 
Dzyaloshinskii-Moriya interaction and the standard su- 
perexchange (the latter active only when multiferroic ma- 
terials that can be controlled by electric fields are in- 
volved) could induce the exchange bias phenomenon at 
FM/G-AFM perovskite oxides interfaces, even when the 
antiferromagnetic spins are compensated. The common 
precondition for the existence of these two mechanisms 
is the presence of oxygen octahedral tiltings at the in- 
terface. Our model highlights the interactions between 
magnetism and lattice distortions, and proposes mecha- 
nisms to understand the exchange bias in FM/G-AFM 
oxides heterostructures. 
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